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Abstract 

This paper reports the analysis of two Chandra X-ray observations of the young 
magnetic B star p Ophiuchus Si. X-ray emission from the star was detected in 
both observations. The average flux is almost the same in both, but during each 
observation the flux shows significant time variations by a factor of two on timescales 
of 20-40 ksec. Each spectrum can be fit by either an absorbed power law model with 
a photon index of ~ —3 or a thin-thermal plasma model with a temperature of ~2 
keV and an extremely low metal abundance (< 0.1 solar). The spectrum of the first 
observation has an apparent line feature at about 6.8 keV, which likely corresponds 
to highly ionized iron Kcr. In contrast, the spectrum of the second observation shows 
an anomalous edge absorption component at E ~-l keV. The continuum emission 
and log (Lx/ Li, 0 i)~ —6 are similar to those of young intermediate-mass stars (Herbig 
Ae/Be stars) although the presence of the magnetic field inferred from the detection 
of non-thermal radio emission has drawn an analogy between p Ophiuchus Si and 
magnetic chemically peculiar (MCP) stars. If the X-ray emission is thermal, the 
highest plasma temperature observed is too high to be explained by the conventional 
theories of magnetic stars, and favors some kind of magnetic dynamo activity, while 
if the emission is nonthermal, it might be related to mass infall. The 6.8 keV line 
and 4 keV edge features are marginal but they give important information near the 
stellar body if they are real. Their physical interpretation is discussed. 

Key words: stars: individual (p Ophiuchus Si), stars: abundances, stars: chem- 
ically peculiar, stars: magnetic fields, X-rays: stars 
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which probably comes from gyro-synchrotron particles in a large magnetosphere, suggests that 
it is a magnetic Bp star. Si has many characteristics of youth; a class III object with a double- 
peaked spectral energy distribution (SED) (e.g. Ward-Thompson 1993, Wilking et al. 2001), 
possession of a compact HII region (A88) and proximity to a plausible star forming cloud SMI 
(Motte et al. 1998). Si may have dissipated most of its disk (< 2.3 x 1O _3 M 0 , Niirnberger et 
al. 1998), so that it should be near or already on MS. 

Observations of Si with Einstein, ROSAT and ASCA have shown relatively strong X- 
ray emission for B stars (Lx ~ lO 30-31 ergs s -1 , Montmerle et al. 1983, Casanova et al. 1995, 
'Kamata et al. 1997), but those observations did not derive timing and spectral properties 
due to the limited photon statistics and severe contamination by a nearby source. This paper 
reports the first X-ray properties for the star, derived from the Chandra data and allows us to 
constrain the emission mechanism for the first time. 

2. Observations and Data Reduction 

Si was observed twice with the Chandra X-ray observatory in timed event mode with the 
onboard CCD camera ACIS (Weisskopf et al. 1996). The first observation (Obsl) was made 
on April 13, 2000 for about 100 ksec. The telescope optical axis on the AC1S-I array pointed at 
the p Ophiuchus cloud core F (c* 2 ooo = 16 h 27 m 18?l, <5 20 oo = — 24°34'21/9: J2000). Si was 14/8 
off-axis and was placed on the ACIS-S3 chip. The second observation (Obs2) was made on May 
15, 2000 for about 96 ksec, pointing at the p Ophiuchus A cloud ( 0:2000 = 16 h 26 m 35?3, £2000 = 
— 24°23'12."9). Si was 0/3 off-axis and was placed on the ACIS-I3 chip. We utilize the level 
2 screened event data, which were processed at the Chandra X-ray Center (CXC) (processing 
software, ver. R4CU5UPD13.2 for Obsl, ver. R4CU5UPD13.2 for Obs2). The reduction and 
further analyses are performed with the software packages CIAO 2.1.3 and FTOOLS 4.2. 

3. Analyses and Results 

3.1. Source Detection & Event Extraction 

In both observations, a bright X-ray source was detected at the optical position of SI 
(error circle, Obsl: ~ 15", Obs2: ~ 0.5") by using the wavdetect package (figure 1). In Obs2, a 
circle of radius of ~ 1.3" includes 95% of the source photons, while in Obsl the 95% radius is 
~ 35" because of the large off-axis angle. In the Obs2 image, no other X-ray source is detected 
within the psf of Obsl. The X-rays in both observations should therefore come from the same 
source uncontaminated by nearby sources. The coordinates derived with the satellite attitude 
data have a systematic offset 1 , which we corrected by cross-correlation of the Chandra detected 
sources with 2MASS counterparts [(Aa, AS) = (0.0", 1.0")]. After the correction, the position 
of the X-ray source is ( 02000 1 ^ 2000 ) = (16 /l 26 TO 34.21 s , -24 d 23'28.2"). The 2MASS position of Si 


1 http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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an absorbed IT model at 96% confidence level (x^/d.o.f = 129.3/102 for a thermal model) 
due to a large residual above 4 keV (table 3). The residual suggests an edge at ~4 keV and a 
gradual rise up to 10 keV. A thermal model including an edge feature reduces the x 2 value to an 
acceptable range (x 2 /d.o./ = 116.5/100 for a thermal model; figure 5). The spectral parameters 
except for the plasma temperature are almost as same as those in Obsl. The column density 
Nh is ~2xl0 22 cm -2 , that is consistent with Ay of Si (~11.7) assuming the Ah— Ay relation 
in the p Ophiuchus cloud (Imanishi et al. 2001). The metal abundance is quite low so that the 
spectrum can be fit by an absorbed power law model with an edge component in the same way. 
As far as we are aware, such an edge feature has not been seen in any stellar X-ray spectra. We 
thus investigate other possibilities for making this feature: a multi-temperature thermal model 
and a single- temperature thermal model with an emission feature near 3.65 keV. However, a 
two temperature model without an edge docs not reproduce the shape of this dip. A gaussian 
line added at 3.65 keV in substitution for the edge feature does not reproduce the excess around 
3 keV. The data do not suffer severe event pile-up nor do the surrounding sources show any 
similar feature. Though our data processed by CXC do not completely recover the effects of 
charge transfer inefficiency(CTI), these effects could not produce the edge feature because they 
continuously affect the pulse height distribution of spectra. There is no plausible instrumental 
origin for this feature to our knowledge. 


4. Discussion 

4-1. The General Characteristics of the X-ray Emission 

The X-ray properties of Si, relatively high plasma temperature and X-ray variability, 
are similar to low-mass young stars rather than early type MS stars. This makes us suspect 
that the X-ray emission might not come from Si itself but from a T-Tauri companion star. 
In fact, Si has a faint close companion (A' = 8. m ) at a projected separation of ~ 0.02", whose 
spectral type is unknown (Simon et al. 1995). However, more than 90% of optically selected 
low-mass stars have X-ray luminosities less than 10 3 ° ergs s -1 (Neuhauser et al. 1995). It is 
therefore likely that the X-ray emission comes from Si itself. 

We here compare the X-ray characteristics of Si with MCP stars and HAeBe stars (or 
young MS stars). Si has —6.5 < log (Lx/Lboi) < —5.5 in the ROSAT band. The soft X-ray 
luminosity is larger than the typical value of He-rich stars (\ogLx / Lboi ~ —7) that are Bp stars 
with strong magnetic fields, and is closer to the X-ray luminosities of non-magnetic Bp stars or 
Ap stars (which have log Lx/Lboi < —6, Drake et al. 1994). In contrast, the X-ray luminosity 
of Si is within the range of HAeBe stars (which have log Lx/Lboi < —4, Zinnecker & Preibisch 
1994; Hamaguchi 2001). The plasma temperature of Si (~2keV) is larger than that of MCP 
stars (<lkeV) but it is typical of temperatures of HAeBe and young MS stars (Hamaguchi 
2001, Feigelson et al. 2002). MCP stars were observed by ROSAT , which was less sensitive 
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4-2.2. Non- Thermal Emission 

Synchrotron radiation produces nonthermal X-rays. It requires at least ~10GeV elec- 
trons to emit X-rays in a magnetic field of a few hundred Gauss, but normal stars are not 
generally thought to produce such electrons. The mechanism therefore has not been considered 
effective in producing non-thermal X-ray emission in stars. Si has nonthermal radio emission, 
which comes from the gyro-synchrotron electrons (A88). We then extrapolate the radio spec- 
trum to the X-ray band. The radio spectrum has a convex shape, which they interpret as the 
free-free absorption in the lower frequency by the thermal X-ray plasma. The linear extrapola- 
tion to high frequencies of the best-fit model of A88 figure 4 is well below the observed X-ray 
emission level (figure 6). Thus the X-ray emission is not explained by the gyro-synchrotron 
emission. 

The radio observation suggests the presence of MeV electrons around Si (ASS). High 
energy electrons which hit a dense region such as the stellar surface, could produce observable 
nonthermal bremsstrahlung X-rays. A similar mechanism has been seen in solar flares. The 
infalling matter is accelerated by magnetic reconnection and collides with the stellar surface 
(see Sakao et al. 1998). The thermal component in solar flares is larger than the nonthermal 
component below 20 keV, but in the case of Si, the nonthermal component should be dominant 
over the thermal component above ~1 keV. The produced plasma could be cooler than the 
thermal plasma of the sun by some reason, for example slower infall speed. The soft thermal 
X-ray emission may be hidden by circumstellar absorption (see subsection 4.3.2). 

4-3. Specific Features Seen in the Spectra 

4-3.1. Line Feature near 6.8 keV 

The apparent line feature seen in Obsl is reproduced with helium or hydrogen like iron 
equivalent to ~lsolar abundance (Ah ~3.1x10 22 cm -2 , kT ~1.2 keV). This model requires 
that only Fe is near solar abundance and the other elements like Mg, Si and S are less than 
solar abundance. Mg and Si have almost the same first ionization potential (FIP) as Fe, and 
therefore the abundances of Mg, Si and Fe should be similar if FIP or inverse-FIP effects work 
(e.g. Gudel et al. 2001). Other types of selective elemental enhancement mechanisms may be 
required. 

The excess is not prominent in Obs2. This could be due to a time variation but it might 
be due to the edge absorption (see the next subsection). Assuming a line at 6.77 keV in Obs2, 
the flux upper limit, 1.6 x 10 -6 photons cm -2 s -1 (equivalent width ~1.1 keV), is consistent 
with the line flux seen in Obsl. 

4-3.2. Edge Feature 

The threshold energy of the edge (3.84 keV < E e d ge < 4.09 keV) includes the K shell 
binding energies of Cl, Ar, K and Ca in neutral or ionizing states (Lotz 1968). Ar has the 
highest relative abundance among them. The abundance of Ca is one third that of Ar, and Cl 
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a nonthermal origin (see subsection 4.2.2). The emitting region and absorber might be on the 
magnetic poles. Assuming the stellar surface density for n a (, s (~10 12 cm -3 ), the geometrical 
depth of the absorber (l a bs) is ~10 n cm (~l/3 i?„). The absorber might be elongated along 
the magnetic field lines. 

5. Probable Models of the X-ray Emitting Region 

We here speculate about the X-ray emission mechanism and the structure around the 
X-ray emitting region. The 6.8 keV line feature seen in Obsl favors thermal emission, while 
the edge feature seen in 0bs2 favors nonthermal emission if it originates from Ar. The spectral 
shape does not change between Obsl and 0bs2 except for these features, and therefore it is 
less conceivable that the emission mechanism changes between the observations. If the X-ray 
emission is thermal, the plasma might be produced by mass infall on the stellar surface. In this 
case, the edge feature should originate from Ca. An emitting plasma with low metal abundances 
and a Ca rich cold absorber could be produced by some elemental selection mechanism. If the 
edge is produced by Ar, this suggests that the X-ray emission is nonthermah However, in this 
case the 6.8 keV line feature seen in Obsl is hard to explain if this feature is real. Nonthermal 
bremsstrahlung emission might be produced by the collision between the surface and infalling 
material, with free fall velocity (~700 km s -1 ). As the heated gases evaporate, they would fill 
the bottom of the polar cap and might be seen as the warm absorber of kT ~0.4 keV. The edge 
absorption might be observed during Obs2 if the warm plasma is moved into the line of sight 
by the rotation of an oblique magnetic field. 

6. Conclusion 

The magnetic B star Si is identified as an X-ray source with large X-ray luminosity (log 
Lx [ergs s -1 ] ~ 30.3) with a precision of ~0"5. The X-ray spectrum implies a large extinction 
of Nh ~10 22 cm -2 , that is consistent with A v of Si. This observation gives good supporting 
evidence that SI itself, i.e. an intermediate-mass star, emit X-rays. The X-rays do not show 
the characteristics of normal early type MS stars; somewhat larger Lx/Lboi ratio (~ —6), small 
but significant X-ray time variations and significantly high plasma temperature (~2 keV). The 
spectra do not show emission line features (except for an anomalous excess at ~6.8 keV in 
one spectrum), while one spectrum shows an edge like feature at around 4 keV. The lack of 
apparent line features makes the limit on the plasma abundance less than 0.1 solar, that is 
quite small compared with other stellar X-ray spectra usually showing ~0.3 solar. 

The plasma temperature and X-ray luminosity are similar to HAeBe and young MS OB 
stars rather than MCP stars. The X-ray emission of Si might be related to its youth. The 
magnetic B star could have an X-ray emission mechanism like magnetic activity or mass infall. 
The small abundances may require a mechanism for selection of elements related to magnetic 
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Table 1. Fitting results of the light curves 


Observations 



1st 

2nd 

Model 


Cons. 

Cons. -J- Exp. 

Cons. 

Mean 

[10 2 cnts s x ] 

2.1 

1.4 (0. 5-1.8) 

2.6 

e-folding time 

[10 4 s] 


3.8 (2.2-8. 3) 

... 

* 2 /d.o.f. 

7-t- ry ' , — 

— r — i — r- : 

154.9/51 

-■ — I ii 

49.7/49 

95.9/49 


Cons.: Constant model, Exp.: Exponential model 

The errors listed in parenthesis quoted for 90% confidence. 


Table 2. Fitting results of the spectrum in the 1st observation (Obsl) 


Model 


IT 

IT -f gauss 

power 

N H 

[10 22 cm- 2 ] 

1.9 (1.7-2. 1) 

1.9 (1.7-2. 1) 

2.4 (2. 1-2.7) 

kT/T 

[keV]/ 

1.7 (1.4- 1.9) 

1.6 (1.4-1. 9) 

3.4 (3. 1-3.6) 

Abundance 

[solar] 

0.0 (0.0-0.11) 

0.0 (0.0-0.09) 


Emission measure 

[10 53 cm" 3 ] 

2.8 (2.2-3.4) 

2.8 (2. 2-3.5) 


Line center 

[keV] 


6.77 (6.56-6.99)* 


Line flux 

[10 -6 photon cm -2 s -1 ] 


3.1 (1. 1-5.0) 


L x (0.5-10 keV)** 

[10 3 ° ergs s -1 ] 

1.8 

1.9 

4.8 

X 2 / d.o.f 


118.5/106 

111.5/104 

116.5/107 


Line width [a) of the Gaussian component is fixed on zero. 

** Other parameters except for the line normalization are temporary frozen in the error estimate. 
* * * Absorption corrected X-ray luminosity assuming the distance of 120 pc. 


Table 3. Fitting results of the spectrum in the 2nd observation (Obs2) 


Model 


IT 

IT x edge 

power 

power x edge 

N . H 

[10 22 cm" 2 ] 

2.0 

1.8 (1. 6-2.0) 

2.5 

2.2 (2. 1-2.4) 

kT/r 

[keV]/ 

1.9 

2.5 (2.1-3. 1) 

3.2 

2.7 (2.5-3. 0) 

Abundance 

[solar] 

0.21 

0.14 (0.0-0.28) 

... 

... 

Emission measure 

[10 53 cm- 3 ] 


1.6 (1.3-1. 9) 



Threshold energy* 

[keV] 


3.96 (3.84-4.07) 


4.00 (3.89-4.09) 

Absorption depth** 



0.53 (0.28-0.81) 


0.67 (0.48-0.90) 

L x (0.5-10 keV)*** 

[10 3 ° ergs s -1 ] 


1.5 

... 

2.6 

X 2 / d.o.f 


129.3/102 

116.5/100 

132.3/103 

111.3/101 


• Threshold energy of the edge component. 

** Absorption depth at the threshold of the edge component. 

• * • Absorption corrected X-ray luminosity assuming the distance of 120 pc. 
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1st Obs. : Mekal * wabs 
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2nd Obs. : Mekal * wabs 



channel energy (keV) 





2nd Obs. : Mekal * edge * wabs 



channel energy (keV) 
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